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We present a search for Wbb production in pp collisions at y/s = 1.96 TeV in events containing one 
electron, an imbalance in transverse momentum, and two 6-tagged jets. Using 174 pb -1 of integrated 
luminosity accumulated by the D0 experiment at the Fermilab Tevatron collider, and the standard- 
model description of such events, we set a 95% C.L. upper limit on Wbb production of 6.6 pb for b 
quarks with transverse momenta p h T > 20 GeV and bb separation in pseudorapidity— azimuth space 
AlZbb > 0.75. Restricting the search to optimized bb mass intervals provides upper limits on W H 
production of 9.0—12.2 pb, for Higgs-boson masses of 105—135 GeV. 

PACS numbers: 13.85Qk,13.85.Rm 



The Higgs boson is the only scalar elementary particle 
expected in the standard model (SM). Its discovery would 
be a major success for the SM and would provide further 
insights into the electroweak symmetry breaking mecha- 
nism. The constraints from precision measurements |l| 



favor a Higgs boson sufficiently light to be accessible at 
the Fermilab Tevatron collider. Although the expected 
luminosity necessary for its discovery is higher than ob- 
tained thus far, the special role of the Higgs boson in the 
SM justifies extensive searches for a Higgs-like particle in- 



4 



dependent of expected sensitivity. Such studies also pro- 
vide an opportunity to investigate the main backgrounds, 
and in particular the interesting and thus far unobserved 
Wbb production process. 

In this Letter, we present a search for a Higgs (H) 
boson with mass m# between 105 and 135 GeV, in the 
production channel pp — > WH — > evbb, at i/s = 1.96 
TeV. The expected WH cross section is of the order of 
0.2 pb for this mass range 0. Our search is based on an 
integrated luminosity of 174 ± 11 pb -1 accumulated by 
the D0 experiment during 2002 and 2003. 

The experimental signature of WH — > evbb relies on a 
final state with one high px electron, two b jets and large 
imbalance in transverse momentum (Et) resulting from 
the undetected neutrino. The dominant backgrounds to 
WH production are from Wbb, ti and single top-quark 
production. The signal to background ratio is improved 
by requiring exactly two jets in the final state, because 
the fraction of ti events that contain at most two recon- 
structed jets is small. We use the high statistics W+ > 2 
jets data to check the validity of our simulation, but re- 
strict the selection to W + 2 6 jets for the final results. 

The D0 detector includes a magnetic tracking sys- 
tem surrounded by a uranium/liquid- argon calorimeter, 
which is enclosed in a muon spectrometer. The tracking 
system consists of a silicon microstrip tracker (SMT) and 
a central fiber tracker (CFT), both located within a 2 T 
superconducting solenoidal magnet Q. The SMT and 
CFT have designs optimized for tracking and vertexing 
capabilities for pseudorapidities < 3 and < 2.5, re- 
spectively . The calorimeter has a central section (CC) 
covering rj up to \r}\ « 1.1, and two end calorimeters (EC) 
extendingcoverage to \rj\ w 4.2, each housed in a separate 
cryostat |5j. For particle identification, the calorimeter 
is divided into an electromagnetic (EM) section, followed 
by fine (FH) and coarse (CH) hadronic sections. Scin- 
tillators between the CC and EC cryostats provide addi- 
tional sampling of developing showers for 1.1 < |ry| < 1.4. 
The muon system consists of a layer of tracking detectors 
and scintillation trigger counters in front of 1.8 T toroids, 
followed by two similar layers behind the toroids, which 
provide muon tracking for |?7| < 2. The luminosity is 
measured using scintillator arrays located in front of the 
EC cryostats, covering 2.7 < |r?| < 4.4. 

Event selection starts with the requirement of an iso- 
lated electron, with px > 20 GeV, in the central region of 
\rj\ < 1.1, but away from boundaries of calorimeter mod- 
ules at periodic azimuthal angle (ip) values @ . Such elec- 
trons are required to trigger the event. The average trig- 
ger efficiency is (94 ±3)% for W + 2 jets events. Electron 
candidates are selected by requiring: (i) at least 90% of 
the energy in a cone of radius ATI = y/(Aip) 2 + (Arj) 2 = 
0.2, relative to the shower axis, is deposited in the EM 
layers of the calorimeter, i.e., EM fraction emf> 0.9; (ii) 
isolation, i.e., that the total energy in a cone of A1Z < 0.4 
centered on the same axis does not exceed the recon- 



structed electron energy by more than 10%; (iii) that the 
energy cluster of the electron candidate has the charac- 
teristics of an EM shower, as determined by the stan- 
dard D0 shower-shape criteria 6]; (iv) that there is a 
track pointing to the EM cluster. These four criteria de- 
fine the initial electron candidates. Electron selection is 
further refined using an electron likelihood discriminant 
based on the above estimators, as well as on additional 
tracking information. The combined reconstruction and 
identification efficiency is determined from & Z —* e + e~ 
sample to be (74 ± 4)% per electron. 

To select W bosons, we require Et > 25 GeV. Events 
with a second isolated lepton (e or (i 7]) with pt > 15 
GeV and \rj\ < 2.4 are rejected to suppress Z + jets and ti 
backgrounds. Only events with a primary vertex at \z\ < 
60 cm relative to the center of the detector are retained. 
At least two jets with px > 20 GeV and \r]\ < 2.5 are 
then required. A jet is defined as a cluster of calorimeter 
towers within a radius A1Z — 0.5 8], having: (i) 0.05 < 
emf < 0.95; (ii) less than 40% of its energy in the CH 
section of the calorimeter; (iii) a distance A1Z to any ini- 
tial electron candidate greater than 0.5. The average jet 
reconstruction and identification efficiency is (95 ± 5)%, 
as determined from 7+ jet events. For selecting b jets, we 
use an impact-parameter based algorithm which has 
been cross-checked with a secondary-vertex reconstruc- 
tion algorithm. 

To improve calorimeter performance, before recon- 
structing the calorimeter objects, we use an algorithm 
that suppresses cells with negative energy (originating 
from fluctuations in noise) and cells with energies four 
standard deviations below the average electronics noise 
(<7„), when they do not neighbor a cell of higher energy, 
E > 4er„. The EM scale is calibrated using the peak 
in the Z — > e + e~ reconstructed mass, and jet energies 
are then corrected to the EM scale using 7+jet events. 
These energy corrections, and the transverse momenta 
of any muons in the event, are propagated into the cal- 
culation of the ^r, which is estimated initially using all 
(unsuppressed) calorimeter cells. 

The D0 detector simulation based on GEANT |J3 an d 
the reconstruction and analysis chain used for data are 
also used for obtaining expectations from the standard 
model, which are normalized to cross sections measured 
in data, or to calculations when no such measurements 
are available. Small additional energy smearing in Et 
and in the energy of the simulated electrons is used to 
obtain better agreement between data and simulation. 

Before applying b tagging, we expect to have two main 
components in the data: W+ jets events and multijet 
events in which a jet has been misidentified as an elec- 
tron (called QCD background in the following). W+ 
jets events are simulated u sing the leading-order matrix- 
element program ALPGEN [ill for the Wjj process (i.e. 
production of W + 2 partons, which are in our case glu- 
ons or u,d,s,c quarks, since the Wbb is simulated sepa- 
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FIG. 1: Distribution of the dijet invariant mass of W + 2 
jets events, compared with cumulative contributions from 
the QCD background (derived from data), the simulation 
of W±jets events and the other SM backgrounds, which 
are small before b tagging. Uncertainties on the simulation 
from systematics of the jet energy scale are indicated by the 
hatched bands. The simulated contributions are normalized 
to the integrated luminosity of the data. 

rately), followed by PYTHIA [l2j for parton showering and 
hadronization. The QCD background is estimated from 
data using measured probabilities for jets to be misiden- 
tified and accepted as electrons. 

The distribution of the dijet invariant mass in W + 2 
jets events is shown in Fig. ^ where it is compared to 
expectation. The Wjj expectation is normalized to the 
data using the next-to-leading-order (NLO) MCFM calcu- 
lation 13], providing a simulated rate for W + 2 jets 
events in agreement with the measured rate. Taking 
into account uncertainties originating from the jet energy 
scale, the shape of the distribution is also well described. 
The systematic uncertainty associated with the selection 
of exactly two jets in the final state has been studied in 
data and in simulations. The rates for W + 3 jets and 
W + 4 jets events, after normalizing to the W + 2 jets 
sample, are described by the ALPGEN and pythia simu- 
lation to within 15% and 6%, respectively. The resulting 
systematic uncertainty on the expectation is ±5%. 

To search for Wbb final states and to suppress back- 
ground, we apply the ^-tagging algorithm to jets having 
at least two tracks, with ^ ackli ^ > 1.0(0.5) GeV. These 
requirements have a typical efficiency per jet of 80% for 
multijets events, which is reproduced to within 5% by 
the simulation. The ^-tagging algorithm uses a lifetime 
probability that is estimated from the tracks associated 
with a given jet. A small probability corresponds to jets 
having tracks with large impact parameters that charac- 



terize 6-hadron decays. Requiring a probability smaller 
than 0.7%, yields a mistag (tagging of u, d, s or gluon 
jets as b jets) rate of (0.50±0.05)%. The tagging effi- 
ciency for a central b jet with p^ between 35 and 55 GeV 
is measured to be (48±3)%. 

The tagging efficiency in the simulation is adjusted to 
the one measured in data. A study of the pr and rj de- 
pendence in data and in simulation indicates a system- 
atic uncertainty on tagging efficiencies of ±6%. When 
tagging light quarks, there is a larger systematic uncer- 
tainty on the efficiency (±25%) that originates from the 
direct application of the algorithm to simulated events. 
This has only a small effect on the final results, since the 
fraction of events with two mistagged jets is < 10% of 
the total number of W + 2 6-tagged jets. For the tagging 
efficiency of c quarks, we use the same data/simulation 
efficiency ratio as for b quarks. 

To reduce the presence of b jets from gluon splitting, 
and to help assure an unambiguous determination of jet 
flavors in simulation, we require the separation between 
the two reconstructed jets (All) to be greater than 0.75. 
In Fig. [3 we show the distribution of the dijet mass for 
W+2 jets events in which at least one jet is b tagged. The 
data are well described by the sum of the multijet back- 
ground and simulated SM processes (cf. Table [!}. The tt 
contribution is simulated with pythia (a t i = 6.77± 0.42 
pb 01)- Single-top production (aw*^tb = 1-98 ± 0.32 
pb, (T g w-^tb = 0.88±0.13 pb 0) is generated with com- 
phep [lfj, assuming a top-quark mass of 175 GeV, and 
is shown in Fig. in combination with other processes: 
Z — > ee, W — ► tv and WZ(—> bb), which are simu- 
lated using pythia with cross sections of 255 pb [l7j . 
2775 pb [21, and 0.6_pb [3, respectively. As for the 
Wjj process, the Wbb contribution is simulated using 
ALPGEN and PYTHIA, requiring p h T > 8 GeV and AlZbb > 
0.4 at the parton level, with o wh i =3.35 pb computed 
at NLO using the MCFM program. WH production is 
simulated with pythia using the computed cross section 
at NLO, which depends on tuh Q- 

To further improve signal/background, we select events 
in which a second jet is b tagged. The final results for 
the number of observed and expected events are given in 
Table ID Data fr om the last column are not used in the 
analysis, but provide a check of the accuracy of our ex- 
pectations for events with two 6-tagged jets in the control 
sample of W+ > 3 jets events, which is dominated by tt 
production. 

The distribution of the dijet mass for events with two 
6-tagged jets is shown in Fig. EH The expected number 
of events is 4.4 ±1.2, of which 1.7 events are expected 
from Wbb production. The dominant systematic uncer- 
tainties on the expectation come from uncertainties on 
the 6-tagging efficiency (11%) and jet energy corrections. 
The uncertainty on the latter propagates to uncertain- 
ties of 7% on Wbb production, 4% on single-top and WH 
production, and 3% on tt production. The total system- 
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atic uncertainty on the expectation is 26%, including the 
uncertainties on cross sections and luminosity (18% and 
6.5%, respectively). 



- D0 W +> 1 b-tagged jets 
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□ W + jets 
I □QCD 
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DijetMass (GeV) 

FIG. 2: Distribution of the dijet invariant mass for W + 2 
jets events, when at least one jet is b tagged, compared to 
expectation (cumulative). The other SM backgrounds include 
single-top events. The simulated contributions are normalized 
to the integrated luminosity of the data. 



section, for p b T > 20 GeV and ATZ bb > 0.75 18]. The 
limits on the cross sections are obtained using a Bayesian 
approach that takes account of both statistical and 
systematic uncertainties. 

The expected contribution from the bb decay of a SM 
Higgs boson, with m H = 115 GeV produced with a W, 
is also shown in Fig. [3 and amounts to 0.06 events. The 
mean and width of a Gaussian fit to this expected contri- 
bution in the mass window 85-135 GeV are 110 and 16 
GeV, a relative resolution of (14 ± 1)%. Similar resolu- 
tions are obtained for Higgs-boson masses in the 105-135 
GeV region. 

No events are observed in the dijet mass window of 85- 
135 GeV. The expected SM background (including Wbb) 
is 1.07 ± 0.26 events, and the expected WH signal is 
0.049 ± 0.012 events, with a signal efficiency of (0.21 ± 
0.03)%. In the absence of a signal, we set a limit on the 
cross section for a(pp — ► WH) x BR(H — > bb) of 9.0 pb 
at the 95% C.L., for a 115 GeV Higgs boson. 

The same study was performed for m# = 105, 125 and 
135 GeV, for which 0, and 1 event were observed in the 
corresponding mass windows. The resulting limits (11.0, 
9.1 and 12.2 pb, respectively) are compared to the SM 
expectation in Fig. ^ and to the results published by the 
CDF collaboration, using a smaller integrated luminosity 
of 109 pb -1 at y/s = 1.8 TeV, but for combined e and /i 
channels |2(ij . 
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FIG. 3: Distribution of the dijet invariant mass for W + 2 b- 
tagged events, compared to expectation (cumulative). The 
simulated contributions are normalized to the integrated lu- 
minosity of the data. The expectation for a 115 GeV Higgs 
boson from WH production is also shown. 

Assuming that the six observed events are consistent 
with the SM, without contributions from Wbb and WH, 
and using the Wbb signal efficiency of (0.90 ± 0.14)%, we 
set a 95% C.L. upper limit of 6.6 pb on the Wbb cross 




• D0, 174pb , WH^evbb 
° CDF Run 1, 109 pb" 1 , e+n 



^daitf mode/ 



105 110 115 120 125 130 135 

c Higgs Mass (GeV) 

FIG. 4: 95% C.L. upper limit on a{pp -> WH) x BR(H -» bb) 
compared to the SM expectation at y/s — 1.96 TeV, and to 
CDF results |H, which were obtained at y/s = 1.8 TeV. The 
predicted WH cross section at 1.96 TeV is approximately 15% 
larger than that at 1.8 TeV. 

In conclusion, we have performed a search for the Wbb 
final state, and have set an upper limit of 6.6 pb on this 
largest expected background to WH associated produc- 
tion. We have studied the dijet mass spectrum of two 
6-tagged jets in the region where we have the best sensi- 
tivity to a SM Higgs boson, and for Higgs-boson masses 
between 105 and 135 GeV we set 95% C.L. upper lim- 
its between 9.0 and 12.2 pb on the cross section for WH 
production multiplied by the branching ratio for H — > bb. 
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TABLE I: Summary for the e+^r+jets final state: the numbers of expected W+ > 2 jets and W + 2 jets events, before 
and after b tagging, originating from WH (for mu = 115 GeV), WZ, Wbb, top production (ti and single-top), QCD multijet 
background, and W or Z+jets (excluding Wbb which is counted separately) are compared to the numbers of observed events. 
The last column shows the same comparison for the control sample of W+ > 3 jets events that contain two 6-tagged jets. 





W+ > 2 jets 


W + 2 jets 


W+ 2 jets 
(1 6-tagged je 


W + 2 jets 
;) (2 6-tagged jets) 


W+ > 3 jets 
(2 6-tagged jets) 


WH 


0.6 ± 


0.1 


0.4 ± 


0.1 


0.14 


± 


0.03 


0.056 


±0.013 


0.015 


±0.004 


WZ 


1.4 ± 


0.3 


1.2 ± 


0.3 


0.38 


± 


0.09 


0.13 


± 


0.03 


0.02 


± 


0.01 


Wbb 


24.7 ± 


6.2 


21.4 ± 


5.3 


6.6 


± 


1.5 


1.72 


± 


0.41 


0.37 


± 


0.09 


tt 


41.4 ± 


8.7 


8.6 ± 


1.8 


2.7 


± 


0.6 


0.78 


± 


0.19 


4.63 


± 


1.11 


Single-top 


11.6 ± 


2.4 


8.3 ± 


1.7 


2.7 


± 


0.6 


0.47 


± 


0.11 


0.30 


± 


0.07 


QCD multijet 


492 ± 


108 


393 ± 


86 


17.1 


± 


4.3 


0.50 


± 


0.20 


0.92 


± 


0.37 


W or Z+jets 


2008 ± 


502 


1672 ± 


418 


43.0 


± 


12.9 


0.78 


± 


0.22 


0.86 


± 


0.24 


Total expectation 
Observed events 


2580 ± 626 
2540 


2106 ± 513 
2116 


72.6 ± 20.0 
76 


4.44 ± 1.17 
6 


7.12 ± 1.89 
7 
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